The length distribution of vortex strings in U(l) equilibrium scalar field theory 
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We characterize the statistical length distribution of vortex strings by studying the non-perturbative 
thermodynamics of scalar (7(1) field theory in 3D. All parameters characterizing the length distri- 
butions exhibit critical behavior at T c and we measure their associated critical exponents. Moreover 
we show that vortex strings are generally self-seeking at finite temperature and that there exists a 
natural separation between long strings and small loops based on their spatial correlation behavior. 
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The accurate knowledge of the thermodynamic den- 
sity and length distribution of a population of topological 
defects is a fundamental ingredient for scenarios of de- 
fect formation and evolution in the early Universe and in 
the laboratory. Although a phcnomcnological model ex- 
ists, due to Kibble JlJ and Zurek J|, for predicting defect 
densities it does not provide us with a microscopic un- 
derstanding of the dynamics of defect populations, their 
length distribution and correlations. These aspects of 
the defect network at the time of formation will have a 
crucial impact on the subsequent network evolution. 

In the past, several attempts have been made at char- 
acterizing the statistics of defect networks. Analytical 
studies of cosmic string distributions at finite tempera- 
ture were based on the assumption that strings form a 
gas of non-interacting random walks ||. The general- 
ization to realistic interactions is problematic and yields 
only qualitative predictions. Numerical support for the 
free random walk picture was found from Vachaspati- 
Vilenkin type algorithms [ Ij . In these simulations phases 
of an underlying field are thrown down at random on sites 
of a lattice and vortex lines are formed by joining points 
around which the phases wind about the circle following 
the shortest path. This yields a (lattice-type dependent) 
fraction of long string of 67 — 80% and a length distribu- 
tion that is Brownian for long enough strings. However 
these studies correspond to the infinite temperature limit 
of the underlying field theory Q , where it becomes triv- 
ially free, and therefore have no direct predictive power 
at criticality, when defects are formed. 

Recently, developments in computer power made pos- 
sible the first studies of defect network properties directly 
from the nonperturbative thermodynamics of field theo- 
ries. These focused their attention on defect densities. In 
this letter we proceed to characterize in detail the equi- 
librium statistics of a network of vortex strings in scalar 
1/(1) field theory in 3D. We measure several universal 
critical exponents associated with the string length dis- 
tribution and construct a detailed picture of the critical 
behavior of the theory in terms of vortex strings. 

The adoption of a dual perspective, i.e. of the descrip- 
tion of the mechanisms underlying criticality in terms of 



non-linear solutions of the theory was found to be the key 
to understanding the Kosterlitz-Thouless || transition in 
2D, which proceeds by vortex-antivortex pair unbinding. 
Generally, given a dual model, critical behavior can be 
described equivalently in either the fundamental or the 
dual picture, the choice between either point of view be- 
ing made based on computational advantage. However, 
the construction of exact dual theories proves very diffi- 
cult in D > 3, except when supersymmetry is unbroken. 
The measurements of the string network universal crit- 
ical exponents presented below provide us with direct 
constraints on dual field theories to t/(l) A0 4 in 3D 0. 

Although U(l) Xcf> 4 is the simplest field theory contain- 
ing vortex lines the model is of direct relevance to the 
thermodynamic description of liquid 4 i/e and of strong 
type II superconductors (such as most new high-T c mate- 
rials) in the absence of an external field. It also embodies 
applications to cosmology, eg., the phase transition asso- 
ciated with the spontaneous breakdown of Peccei-Quinn 
symmetry and the formation of axionic strings j^] . 

The methods used to drive the system to canonical 
equilibrium are familiar from stochastic field quantiza- 
tion. A more detailed description can be found elsewhere 
H . The equations of motion for the fields <fn are 
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where i,j G {1,2} and i ^ j in Eq. (|l]). Ti(x) is the 
Gaussian noise characterized by 

(Ux)) = o, (r j; (x)r,(x')} = - O, (2) 

where x, x' denote space-time coordinates. Eqs. (|l[|) en- 
sure that the fluctuation-dissipation theorem applies in 
equilibrium, which will result for large times at temper- 
ature T = We use a computational domain of size 
N 3 = 100 3 , lattice spacing Ax = 0.5 and impose periodic 
boundary conditions in space. At criticality all universal 
exponents are independent of Ax. We assume ergodic- 
ity of the evolution after equilibrium was reached. Using 
this fact we construct the string network from the phase 
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of the complex field <J> = 4>\ + i(f>2 , at regular time inter- 
vals. In all measurements presented below a number of 
samples between 400 and 1200 was used. This permits 
us studies of the string network with unprecedented high 
statistics. The inverse critical temperature _8 C , measured 
from field correlators, is j3 c — 1-91 ± 0.01 ||. 

We separate strings as has become usual in cosmology 
into two distinct populations, one of loops, vortex lines 
shorter than ~ TV 2 and another of long strings longer 
than any loop. Note that there is a priori no reason for 
this separation or for the definition of the boundary be- 
tween the two classes. First we deal with the statistics of 
loops. Characteristic loop distributions above and below 
the critical temperature are shown in Fig. El. 
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FIG. 1. The loop length distribution for several j3. An 
arbitrary distribution (solid) is plotted for comparison. 

In order to characterize the loop length distribution 
n(Z) quantitatively we fit it for each temperature to an 
expression of the form: 

(n(l)) = A I-t cxp(-/3a cff 0, (3) 

where the average is over the canonical ensemble. cr c s 
is the T dependent effective string tension and has di- 
mensions of mass per unit length. The temperature at 
which it vanishes defines the so-called Hagedorn temper- 
ature Th, the point at which long string appears in the 
system. We present our estimate of errors, mostly re- 
sulting from performing multiparameter fits, in the final 
measurements of universal quantities only. Statistical er- 
rors are generally very small. 7 is related to the string 
interactions. For Brownian (i.e. non-interacting) strings 
7 = 5/2, whereas 7 > 5/2 or 7 < 5/2 for self-avoiding 
or self-seeking walks, respectively. The loop distributions 
measured in equilibrium are in fact very well fit by Eq. ([?]) , 
with A(/3), 7(/3) and cr e g(/3) determined by y 2 minimiza- 
tion. These three quantities undergo radical behavior 



changes near the T c lending support to a recent claim 
II that the critical temperature for the fields, T c and the 
Hagedorn temperature for the strings, Th coincide in the 
infinite volume limit. 

Fig. ^| shows the temperature dependence of 7 and a e g. 
We observe that as the critical point is approached from 
below cr c ff vanishes. A fit to a c ff in this regime gives 

<7 eff (/3) oc {p-p H y st 

with H = 1.96 ±0.01, v st = 1.50 ±0.01 and where H 
17 ; 1 . Note that j3n does not coincide exactly with /3 C . 
This discrepancy was observed previously in [^J , where it 
was shown that as the computational domain is increased 
the temperature at which long strings first appear /3h 
migrates to C . In addition the relative precision of the 
measurement of cr e ff close to /3 C is poor since it requires a 
large number of long loops, lying on the boundary of the 
long string/loop separation, eg., for (i = 1.975, we have 
that I ~ {Po-csT 1 = 0.44 x 10 4 ~ N 2 . Below a more 
accurate measurement of (3h will be obtained from the 
behavior of the length distribution of long strings. 
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FIG. 2. a) ffrffOS) and b) 7 (/3) 

Also very interesting is the value of v> st . This universal 
exponent for <7 e ff in the global U{1) theory, corresponds 
to the exponent t'duai = v s %/2, associated with the corre- 
lation length of its dual field theory H. We believe this 
measurement provides a first direct determination lead- 
ing to the classification of the universality class to which 
the dual formulation belongs. Note that v& ua \ = 0.75 dif- 
fers clearly from the mean-field prediction z^mf = 0.5, for 
any O(N), N > 1, scalar field theory. Remarkably t'duai 
coincides with the exact value for the critical exponent of 
an ensemble of interacting polymers in 2D || . 

As for 7, it is especially interesting to note that 7 < 5/2 
generally at any finite temperature, see Fig. ||. The dis- 
crepancy between our observations and the Brownian re- 
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suit cannot be attributed to small computational domain 
effects. The computation of the Brownian distribution in 
a finite volume generates exponential corrections jnj to 
Eq. (||) , which fail to describe the changes in the distribu- 
tions correctly and are moreover [3 independent. Above 
T c , j(/3) ~ 5/2 - 0.0385 (3 and 7 (/3 c ) = 2.23 ± 0.04, Q. 

Fig. g shows the (3 dependence of A and of the average 
walk step size a(/3), defined in terms of A(j3) and 7 (/3) by 



A = 



3/2 



,7-3 



(4) 



The numerical prefactor in Eq. (|]) is obtained from the 
computation of (n(l)) in the free random walk case [pi. 
We observe that the amplitude of the distribution A((3) 
varies drastically with temperature in the vicinity of T c . 
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FIG. 3. a) Fhe amplitude of the loop distribution A(f3) and 
b) the average step length a{0) and the corresponding best 
fits Eqs. (|^) at criticality. 

A(j3) alone determines the /3-dependent probability of 
finding the smallest loops. In the vicinity of the critical 
point A suffers a dramatic drop until (3 C is hit, followed 
by a sudden increase. At /3 = f3 c dA/d(3 is discontinuous. 
The sharp drop in A signifies that a((3) has grown and 
with it the possibility for producing longer strings. 

From general considerations at a second order phase 
transition we can expect that all characteristic lengths 
will present divergences as T — > T c . a((3) is no exception. 
Note however that while £(/3), the correlation length as- 
sociated with the connected field 2-point function, mea- 
sures isotropic correlations of field excitations, a refers 
to correlations along the strings. It is therefore interest- 
ing to measure the critical exponents associated with a. 
Fig. ||b shows the best fits to a(j3) at criticality: 



a(P) -lac 1/(0 -0 c )"i, f3>0 c 
a(0)-lot:l/(p c -0) v *, /?<&, 



(5) 



with zri~ = v~ = 0.5±0.1. These values are typical 
of mean-field predictions and differ from v — 0.671 [6), 
associated with £. 

Above the critical point the behavior of A is given by 
A(f3) ~ 0.092 + 0.006 (3. Note that in addition, given 7 
and above T c (a c s — 0), the total density of loops 



Ploop 



1=4 



I < n(l) >- 



A 



4(7-2) ( 7 _ 2 ) 



(7 > 2) (6) 



is proportional to A in the infinite volume. For (3 = 0, 
when the distribution becomes Brownian the value of A 
does indeed exactly coincide with our measurement of 
Pioop(0) = A(0) = 0.083 in @. 

In order to complete the picture of the role of vortex 
string fluctuations in the thermodynamics of the theory 
we characterize in Fig. |], the length distribution of long 
strings and show the mean distance between two string 
segments R(n) as a function of the number of steps (in 
lattice units) along the string. The distribution of long 
string, is notoriously difficult to measure because the to- 
tal number of long strings in any computational domain 
is typically only a few. Here we could probe 600-1200 in- 
dependent string networks, a large number compared to 
any study in a computational domain of this size. How- 
ever, in such a set of data the number of long strings of a 
given length is never more than one in the total ensem- 
ble of measurements and a direct fit to the data becomes 
impossible. To obviate this problem we adopted a maxi- 
mum likelihood test, see eg. |l^], to fit the data from all 
samples to a distribution of the form Al 11 . In the free 
random walk approximation 7; = 1, as the result of the 
periodicity of the boundary conditions [|l0| . 

Fig.|]a shows the temperature dependence of 7 j. For 
T > T c , 7; is close to the Brownian value 7; = 1 converg- 
ing to it as (3 — > 0, as 7/ ~ 1 + 0.1147 (3. In the small 
band of temperatures below T c for which we identified 
long string, 7; assumes much larger values. 

This sharp increase in 7; is an indication that the dis- 
tribution no longer has the simple form of a power-law. 
In analogy to the form of the loop distribution it marks 
the onset of exponential suppression. Note that it was 
necessary to look at strings classified as long to notice 
this suppression, allowing us a more precise insight into 
the temperature variation of cr e ff, just below criticality. 
This behavior reinforces our previous conjecture that true 
long string is present in the system only above T c which 
would, in the infinite volume limit, coincide with Tj{. 

A complementary test on the interacting nature of the 
strings is provided by R(n), shown in Fig. Generally 
one expects R(n) ~ n Q , where n is the number of steps 
in lattice units. For Brownian walks it is well known that 
a = 0.5. On the other hand a > 0.5 if the interactions are 
repulsive, (a ~ 0.59 for self-avoiding walks) and a ~ 1/3 
if they are attractive at low T ]l3] ]. 

In U(l) scalar field theory two string segments, at dis- 
tance d interact with potential V(d) ~ nx.n^ ln(r/m), 



for d > 



the width of the string. Here fti is the 
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unit vector normal to the i th string segment times the 
topological charge of the vortex. Thus vortices interact 
with a potential which grows with distance. At short 
distances d < mT x the interactions vanish since the su- 
perposition of two vortices (or of a vortex-antivortex) is 
also a solution of the static field equations. 

Small loops of string are the direct equivalent in 3D 
of vortex-antivortex pairs in 2D, since each segment will 
tend on the average to see, normal to the string, another 
string segment with opposite vorticity. This fact permits 
the nucleation of small loops of string below T c . 
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FIG. 4. a) 7{(/3) for the long string length distribution and 
b) R(n) at ft = 1.8 for loops (A) and long strings (O) an d 
for loops at /3 = 2.05 (□). Solid curves show the behavior of 
free random (a = 0.5) and self-seeking (a — 1/3) walks. 

As T increases longer loops can be nucleated. The na- 
ture of the interactions makes it advantageous to create 
loops where the correlations between segments of oppo- 
site vorticity are maximal. Consequently loops behave on 
the average as random walks for a short distance, while 
the interactions are irrelevant, and beyond become self- 
seeking, see Fig. |]b. This distance increases with T. 

As the critical point is approached the total string den- 
sity mounts to a universal critical value ptot(Pc) — 0.2. 
At this point enough string exists in the system that 
the effects of the interactions shield each others out effi- 
ciently. Precisely then a population of long strings can be 
nucleated, interacting with itself through the background 
loop sea with a potential that is almost vanishing. In this 
sense long strings are born quasi-Brownian. This fact 
is demonstrated in Fig. |]b. After a short transient for 
small distances, where they are self-avoiding (otherwise 
they would have become a loop), long strings assume an 
extraordinarily clear random walk behavior, correspond- 
ing to Hausdorf dimension Dh — 2. This establishes the 
difference between the two string populations, since loops 
remain self-seeking (Dh ~ 3) even for T > T c , see Fig. |. 



In conclusion, we have found quantitative evidence 
that the phase transition in 3D scalar U(l) field theory 
proceeds by vortex string nucleation and in particular is 
triggered by long string proliferation, as often suggested 
]14| . We established that a population of quasi-Brownian 
long strings exists above T c only, and that it is different 
from small loops which remain self-seeking. We char- 
acterized the length distribution of both loops and long 
strings, and measured associated universal critical expo- 
nents. The latter provide the first direct constraints on 
field theories dual to the U(l) scalar model in 3D. The 
critical behavior of the length distributions confirms our 
previous claim that the Hagedorn phase transition in the 
strings coincides with the critical point of the theory. 
The Ginzburg temperature, (flc — 2.34), seems to play 
no particular role in the thermodynamics. 
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